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ABSTRACT 

We present the optical narrow line ratios in an SDSS based sample of 3 175 broad Ha selected 
type 1 AGN, and explore their positions in the BPT diagrams as a function of the AGN and 
the host properties. We find the following: 1. The luminosities of all measured narrow lines 
(Ha,H/3,[0 III],[N II], [S II], [O I]) show a Baldwin relation relative to the broad Ha luminos- 
ity LbHa, with slopes in the range of 0.53 —0.72. 2. About 20% of the type 1 AGN reside 
within the 'Composite' and 'SF' regions of the BPT diagrams. These objects also show ex- 
cess narrow Ha and UV luminosities, for their LbHa, consistent with contribution from star 
formation which dominates the narrow lines emission, as expected from their positions in 
the BPT diagrams. 3. The type 1 which reside within the AGN region in the BPT diagrams, 
are offset to lower [S Il]/Ha and [N Il]/Ha luminosity ratios, compared to type 2 AGN. This 
offset is a selection effect, related to the lower AGN/host luminosity selection of the type 2 
AGN selected from the SDSS galaxy sample. 4. The [N Il]/Ha and [N Il]/[S II] ratios in type 
1 AGN increase with the host mass, as expected if the mass-metallicity relation of quiescent 
galaxies holds for the AGN narrow line region. 5. The broad lines optical Fe II is higher for 
a higher [N Il]/Ha, at a fixed L\, Q \ and Eddington ratio L/Leaa- This suggests that the broad 
line region metallicity is also related to the host mass. 6. The fraction of AGN which are LIN- 
ERs increases sharply with decreasing L/L^. This fraction is the same for type 1 and type 2 
AGN. 7. The BPT position is unaffected by the amount of dust extinction of the optical-UV 
continuum, which suggests the extincting dust resides on scales larger than the NLR. 
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1 INTRODUCTION 

The gas located on 1 - 1 000 pc scale from the center of Active 
Galactic Nuclei (AGN) plays a role in several important processes, 
which are not well understood. This gas is the source of AGN fuel, 
and may absorb AGN energy and momentum output, thus poten- 
tially coupling the growth of the bulge with the growth of the cen- 
tral black hole. It is enriched during the life cycle of stars near the 
nucleus, and therefore traces the star formation history. It also re- 
processes the AGN ionization continuum, which originates from a 
few Schwarzschild radii, and thus its emission may allow to con- 
strain the accretion mode in the innermost regions. 

The most prominent optical signature of the circumnuclear 
gas in AGN is its emission lines, which have widths typical of 
the galaxy potential (~ 300 km s~'). These lines are known as 
the narrow emission lines, and the emitting region as the narrow 
line region (NLR). The vast majority of NLR analyses were per- 
formed on type 2 AGN where the central source is obscured, partly 
because the narrow lines are not blended with the broad emission 
lines, which dominate the emission features in unobscured type 1 
AGN. Previous studies of the NLR of type 1 AGN were either lim- 
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ited to the most prominent forbidden lines (e.g. Boroson & Green 
1992, using [O III] A5007), limited to small samples (e.g. Baldwin, 
Phillips & Terlevich 1981, hereafter BPT, Cohen 1983, Ho et al. 
1997b, Rodrfguez-Ardila et al. 2000, Veron-Cetty et al. 2001, Di- 
etrich et al. 2005), or limited to samples of very weak type 1 AGN 
(e.g. Greene & Ho 2007) in which the narrow lines become more 
prominent (Stern & Laor 2012b, hereafter Paper II). 

Here we provide the first study of the NLR properties of a 
large sample of 3 175 type 1 AGN, defined in Stern & Laor (2012a, 
hereafter Paper I) with minor adjustments below. The Tl samples 
spans a black hole mass range of 10 6 < Mgn < 10 9,5 M Q and a 
bolometric luminosity range of 10 42 < L bol < 10 46 erg s" 1 . In con- 
trast with earlier studies of type 2 AGN, here the innermost AGN 
is directly observed. We use the narrow line measurements, com- 
bined with the AGN SED and broad line measurements, to address 
the following questions: 

How complete is the BPT classification of AGN? The BPT 
diagrams (BPT and Veilleux and Osterbrock 1987, hereafter VO) 
compare the ratio of the [O III] to H/3 luminosity (for brevity 
[Om]/H/3), with [N n]A6583/Hce, [S II](AA6716,6731)/Ha, 
and [O I](A6300)/Ha. These line ratios provide a measure of the 
relative strength of the higher energy ionizing photons, and thus dif- 
ferentiate between stellar and AGN excitation. These diagrams are 
widely used to define type 2 AGN samples, using separation lines 
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based on theoretical models (Kewley et al. 2001, hereafter KeOl), 
and based on the observed distribution of star forming galaxies 
(Kauffmann et al. 2003, hereafter Ka03). 

The BPT/VO AGN selection criteria are commonly viewed as 
necessary and sufficient conditions to define AGN. However, AGN 
samples selected by other means show these selection criteria may 
not be necessary. In a hard X-ray selected sample, a unique signa- 
ture for AGN emission, Winter et al. (2010) found that five out of 
60 objects are in the Star Forming (SFs) galaxies regime, i.e. be- 
low the Ka03 line in the [N Il]/Ha panel of the BPT diagrams, and 
five more are between the Ka03 line and the KeOl line, i.e. 'Com- 
posites'. In the Mbh < 10 6 ' 2 M© type 1 sample of Greene & Ho 
(2007), 39% of the objects are SFs or Composites. This fraction 
dropped to 18% when the spectra was taken from a narrower slit 
(Xiao et al. 2011). On the other hand, only 3% of radio loud AGN 
are classified as Composites or SFs (Buttiglione et al. 2010). Using 
the Tl sample, which is selected independently of the narrow line 
properties, we derive the completeness of the BPT-based selection 
criteria, and its dependence on the AGN emission properties. 

How are the properties of the NLR gas related to AGN and 
host properties? In low z type 2 AGN, the value of [N Il]/Ha, 
which follows NLR metallicity, Z NLR , has been found to modestly 
increase with host mass M„ (Groves et al. 2006) and with host ve- 
locity dispersion a % (Annibali et al. 2010). These trends are associ- 
ated with the known Z relation of quiescent galaxies (Lequeux 
et al. 1979, and citations thereafter). The Znlr — M, relation is 
also implied by the fact that most AGN reside in massive galax- 
ies (Ka03) and have Z NLR > Z Q (Storchi-Bergmann et al. 1998, 
Groves et al. 2004, 2006), while the rare low M» AGN have low 
Znlr (Kraemer et al. 1999, Barth et al. 2008; Ludwig et al. 2012). 
However, these samples are dominated by low Lboi AGN, since they 
are based on the detectability of the host galaxy, and therefore are 
limited to a small volume where luminous AGN are rare. 

In high Lb i AGN at high z, anM, - Z relation can be in- 
ferred from the increase of Zblr with L^i (Hamman & Ferland 
1993, 1999, Nagao et al. 2006a), and a likely relation of L bol -M*. 
Though Znlr an d Zblr are related (Shields et al. 2010), there 
seems to be another variable beyond M* which determines Zblr, 
probably related to the accretion rate in Eddington units (L/L^, 
Shemmer & Netzer 2002, Shemmer et al. 2004, Shields et al. 2010). 
Therefore, it is interesting to compare Znlr with Ax>l directly. Most 
narrow line measurements in high Lboi AGN are based on narrow 
line radio galaxies samples (De Breuck et al. 2000, Vernet et al. 
2001, Iwamuro et al. 2003, Nagao et al. 2006b). These studies 
measured UV line ratios, except Iwamuro et al. which measured 
non-BPT optical line ratios. Comparison of NLR properties de- 
rived from different lines can be ambiguous, due to degeneracies 
in the photoionization models (Nagao et al. 2006b). Therefore, the 
dependence of Z NLR and other NLR properties on L\, \ is still an 
open question. In this work we derive indicators of Znlr based on 
the BPT ratios, for a large dynamical range in Lb i- Using the large 
size of the Tl sample, we also decouple the dependence of Znlr 
on Lboi an d on M», and compare Znlr with Zblr- 

Is the ratio of UV to X ray luminosity a measure of the slope 
of the ionizing spectrum? Due to Galactic absorption, the ionizing 
part of the AGN spectrum in the extreme UV is generally unavail- 
able. Laor et al. (1997) showed that the mean 2 keV luminosity L x 
of PG quasars is consistent with an extrapolation of the mean EUV 
slope (Zheng et al. 1997, Telfer et al. 2002). Therefore, the interpo- 
lated slope between L uv and L x , Oox, may be a good estimate of 
the true ionizing slope. Since the BPT diagrams provide an inde- 



pendent constraint on the ionizing slope, we test this hypothesis by 
comparing the BPT ratios with a ox in the Tl sample. 

A related issue concerns the location of the optically thin dust 
found in type 1 AGN samples (Richards et al. 2003, Gaskell et al. 
2004, Paper I), which can harden a ox - If the extincting dust is lo- 
cated within the NLR, the NLR will see a harder spectrum, and the 
BPT ratios will vary with the amount of reddening. If the extincting 
dust originates from larger scales, the NLR will absorb the original 
ionizing spectrum, and the BPT ratios will remain constant. Be- 
low, we constrain the location of the extincting dust using the BPT 
diagrams. 

Is the Seyfert-LINER transition related to other emission prop- 
erties? Kewley et al. (2006, hereafter Ke06) found a bimodality in 
the BPT diagrams between high ionization Seyferts and low ioniza- 
tion nuclear emission line regions (LINERs, Heckman 1980). They 
showed the Seyfert-LINER transition is related to L/L E( jd, as noted 
previously by Ho (2002). This transition has also been claimed to 
be related to the existence of the broad lines, due to the low de- 
tection fraction of broad lines in LINERs (Ho et al. 1997b). We 
address these suggestions based on the Tl sample. 

The paper is organized as follows. In §2. 1 - §2.3 we summa- 
rize the creation of the Tl sample and the measurement of the AGN 
and host properties, detailed in Papers I and II. In §2.4 we describe 
the comparison type 2 sample we use, and account for differences 
in the measuring procedures. In §3 we extend the relative decrease 
with Lboi (the Baldwin effect) found in Paper II for [O III] and Ha, 
to the Hj8, [N II], [S II] and [O I] lines, first used here. We then pro- 
ceed in §4 to measure the BPT ratios of the Tl sample, and their 
dependence on AGN and host characteristics. In §5, we analyze 
objects which occupy a region in the BPT plots which is not pop- 
ulated in type 2 samples. Analysis of LINERs and Composites is 
performed in §6 and §7. In §8 we use the BPT ratios to constrain 
the AGN ionizing spectrum, and the location of the reddening dust. 
We summarize our results in §9. 

Throughout the paper, we assume a FRW cosmology with SI 
= 0.3, A = 0.7 and H Q = 70 km s" 1 Mpc" 1 . 



2 THE DATA 

2.1 The Tl sample selection 

The Tl sample is selected from the 7 th data release of the Sloan 
Digital Sky Survey (SDSS DR7; Abazajian et al. 2009). The SDSS 
obtained imaging of a quarter of the sky in five bands (ugriz; 
Fukugita et al. 1996) to a 95% r band completeness limit of 22.2 
mag. Objects are selected for spectroscopy mainly due to their 
non-stellar colors (Richards et al. 2002), or extended morphol- 
ogy (Strauss et al. 2002). The spectrographs cover the wavelength 
range 3800A~9200A at a resolution of ~ 150 km s _1 , and are 
flux-calibrated by matching the spectra of simultaneously observed 
standard stars to their PSF magnitude (Adelman-McCarfhy et al. 
2008). 

We use SDSS spectra which have 0.005 < z < 0.31 and are 
classified as quasars or galaxies. To ensure a reliable decomposition 
of the broad and narrow components of Ha, we use only spectra 
with S/N > 10 and a sufficient number of good spectral pixels in the 
vicinity of Ha, as detailed in Paper I. These requirements are ful- 
filled by 232 837 of the 1.6 million spectra in DR7, named the par- 
ent sample. The spectra are corrected for foreground dust, using the 
maps of Schlegel et al. (1998) and the extinction law of Cardelli et 
al. (1989). Each spectrum is then fit with three galaxy eigenspectra 
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Object name 


LbHa 


Av 


M, 


Luv 




J000202.95-103037.9 


41.9 


2310 


11.1 


43.8 


-1.50 


J000410.80- 104527.2 


42.6 


1360 


11.1 


44.6 


-1.57 


J00061 1.55+145357.2 


42 A 


3320 


11.2 


44.0 


-1.57 


J000614.36-010847.2 


41.6 


3910 


10.8 


43.2 


-1.54 


J000657.76+152550.0 


41.5 


3020 


10.2 


43.0 


-1.57 



Table 1. The AGN and host characteristics of the Tl sample objects. The 
values of LbHa and Luv are in log erg s _1 , Av is in km s _1 , and M t is 
in log M.0. The last column lists notes for M„, Luv and Lx, separated by 
commas: 'U' indicates an upper limit, 'N' indicates not available, and 'E' 
indicates a large error. The electronic version includes all 3 175 Tl objects. 

representing the host (see §2.2.4 below), and a Lx °c A~ L5 power 
law representing the AGN continuum. The host is subtracted, pro- 
ducing a spectrum free of stellar absorption features, excluding the 
Balmer absorption lines, which are handled at a later stage (see 
§2.2.4). We also subtract a featureless continuum, derived by inter- 
polating the mean continuum level at 6125A-6250A and 6880A- 
7000A. The residual flux at 6250A-6880A (±14, 000 km s" 1 from 
Ha) is then summed, excluding regions ±690 km s from the 
[O I] 116300,6363, [N 11] 116548,6583, [S 11] 116716,6731 and 
Ha narrow emission lines. We find 6986 objects with significant 
residual flux, which is potentially broad Ha emission. 

For the objects with significant residual near Ha, we fit the 
profiles of the broad and narrow Ha, and the [O III] 15007, [O I], 
[N II] and [S II] doublets mentioned above. Narrow lines are fit 
using 4 th -order Gauss-Hermite functions (GHs; van der Marel & 
Franx 1993) and an up to 10 th -order GH is used for the broad Ha 
profile. Further details are given in §2.4 of Paper I and §§2.3-2.4 
of Paper II. The following criteria are applied to the broad Ha fit, 
in order to exclude objects in which the residual flux is not clearly 
BLR emission: the FWHM (Av) of the fit is required to be in the 
range 1 000 - 25000 km s -1 ; the total flux of the fit, and its flux 
density at the line centre, are required to be significant. As [O III] 
and H/3 are used extensively in this paper, we require them to have 
a sufficient number of good pixels in their vicinity for the fit to be 
reliable, as detailed in Paper I. 

Of the 3 243 objects that pass these criteria, we use here 3 175 
objects in which our algorithm achieved reliable narrow line fits 
(see below). Due to the small relative number of objects in which 
the fitting algorithm did not succeed, we do not attempt to improve 
the algorithm further. The broad Ha luminosity (LbHa) an d Av of 
the 3 175 objects of the Tl sample are listed in Table 1. The selec- 
tion effects implied by our selection criteria are detailed in Paper 
I. 

2.2 Narrow line measurements 

The narrow line luminosities of the Tl sample are listed in Table 
2. Below, we address the limitations of our fitting algorithm, which 
deblends the narrow lines from the broad lines and from the stellar 
absorption features. The success of the deblending can be further 
verified with higher S/N spectra, where the transitions between the 
different components are more prominent. Therefore, we corrob- 
orate our results by analyzing mean spectra of different Tl sub- 
groups, which have an effectively higher S/N. 

2.2.1 Bad pixels 

The main source of bad pixels in the SDSS spectra is poor sky sub- 
traction, which degrades the spectrum mainly at 1 > 8000A. There- 



fore, the [S II] and [O I] lines are not measurable in 612 (19%) and 
190 (6%) of the Tl objects, respectively. These objects are marked 
in Table 2, and are disregarded in figures where the line is used. 
Objects in which one of the other lines used in this work has bad 
pixels do not enter the Tl sample (§2.1). 

2.2.2 Upper limits 

Our algorithm can robustly detect the six different narrow lines if 
their mean flux density Fx is 2-3.5 times the local flux density er- 
ror. The exact value depends on how blended a specific line is with 
other spectral features, and is listed in Table 3. Upper limits on 
fluxes of lines with lower Fx are derived by assuming a Gaussian 
profile, with a flux density equal to the minimum Fx required for 
detection and the width fit to the other narrow emission lines. Ob- 
jects with upper limits are noted in Table 2. 

The Tl sample detection fractions of the different lines are 
listed in Table 3. The detection fractions are all > 77%. 

2.2.3 [O lll]-like narrow lines 

As noted in Papers I and II, in 15% of the sample the fit yielded 
FWHM(nHa) > 1.5x FWHM([0 III]). These objects have non- 
or barely-detectable narrow lines near Ha, and there is no clear 
transition between the broad and narrow components of the Balmer 
lines. Therefore, we fit the narrow lines near Ha in these objects 
with a FWHM, 3 rd and 4 th GH parameters equal to those found for 
[O III]. 

An eye-inspection of the narrow H/3 fits yielded another 188 
objects (6%) without a clear NLR/BLR transition, despite hav- 
ing FWHM < 1.5 x FWHM([0 III]). We refit these objects with 
[O IIl]-like profiles, and updated the relevant narrow line fluxes. 
The new fit failed on 68 of the objects (reduced ■/} > 2). Due to 
their relatively small number, we did not attempt to improve the fit, 
and simply removed these 68 objects from the sample. This change 
in the narrow line fluxes of 6% of the Tl sample has a negligible 
effect on the results presented in Papers I and II. 

The narrow Ha, H/3, and [N II] line fluxes are less certain 
in objects fit with an [O IIl]-like profile. Therefore, throughout the 
paper different symbols are used when these measurements are uti- 
lized. These objects are also noted in Table 2. 

2.2.4 Strong stellar Balmer absorption 

We model the stellar absorption features by fitting the first three 
Yip et al. (2004) eigenspectra (ESa) to the SDSS spectra, together 
with a power law for the AGN continuum. Since the Yip et al. ESa 
have emission lines, in ESI we replace the lines with the absorption 
features of the Hao et al. (2005) ESI (detailed in §2.2 of Paper 
II). This step is justified since both ESl's represent an old stellar 
population. In ES2 and ES3, which represent a younger population, 
an emission line free ES is not available, so we simply interpolate 
over the lines. 

Since the absorption lines are significantly wider than the 
emission lines, an interpolation over the emission will not remove 
the entire absorption feature. However, near Ha the interpolation 
is done also over the [N II] lines which flank Ha. Therefore, our 
fit does not account for the entire Ha absorption feature of young 
stars. In Paper II, we found that in the 5% of the Tl objects that 
have L niia < 3 A x L^(host), the L nlia are underestimated due to 
improper subtraction of the stellar absorption. Now, the narrow H/3 
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stellar 



Object name 


H/3 


[O III] 


Ha 


[Nil] 


[S II] 


[0 i] 


Notes 


robust 


absorption 


BPT-[N II] 


BPT-[S II] 


BPT-[0 I] 


J000202.95-103037.9 


40.7 


41.4 


41.5 


41.3 


41.0 


40.2 




+ 


+ 


AGN 


Seyfert 


Seyfert 


J000410.80-104527.2 


41.3 


41.5 


41.9 


41.6 


-1 


-1 


-,-,-,-,N,N 


+ 


+ 


Composite 


SF 


SF 


J00061 1.55+145357.2 


40.2 


40.6 


41.0 


40.7 


40.5 


39.8 


-,-.-.-,-.u 


+ 


+ 


Composite 


SF 


Seyfert 


J000614.36-010847.2 


40.4 


40.9 


41.1 


40.7 


-1 


39.6 


-,-,-,-,N,- 


+ 


+ 


Composite 


SF 


SF 


J000657.76+152550.0 


40.2 


40.8 


40.7 


40.0 


40.3 


39.7 




+ 


+ 


Composite 


Seyfert 


Seyfert 



Table 2. The narrow line measurements of the Tl sample. All luminosities are in log erg s . Notes on the measurements of the six lines are separated by 
commas in column 8, ordered as in the table. A 'U' designates an upper limit, an 'N' designates bad pixels (for [S II] and [O I] only) which are also marked as 
— 1 in the corresponding luminosity. Objects in which the Ha, H/3, or [N II] narrow line measurements are not robust (§2.2.3), or the Balmer lines are affected 
by strong stellar absorption features (§2.2.4), are marked by an 'x' in the respective following columns. Other objects are marked by a '+'. The last three 
columns list the classification of each object in the corresponding BPT panel. The electronic version includes all 3 175 Tl objects. 



Narrow line X (A) min F% / Detection Fraction 



H/3 
[O III] 
[Ol] 

Ha 
[Nil] 
[S II] 
[S II] 



4961 
5007 
6300 
6563 
6583 
6716 
6731 



3 

3.5 
2.5 

3 
3 

2.5 
2 



0.84 
0.99 
0.77 
0.98 
0.92 
0.91 
0.92 



Table 3. The detection fractions of the narrow lines used in the paper. Col. 3 
notes the minimum flux density required for detection, in units of the local 
flux density error. 



emission line is weaker than Ha, and therefore more suspect to 
significant biases due to improper subtraction of the stellar absorp- 
tion features. However, near H/3 the interpolation in ESI and ES2 
is performed only over the narrow Hj8 line, so the wide part of 
the stellar absorption feature is accounted for by our fit. There- 
fore we mark the same objects as in Paper II, i.e. objects with 
inHct < 3^ x L%(hD8t), as objects with potentially underestimated 
LnHa and AiH0 ■ We verify below this suffices in order to identify 
objects with offset L nll p values. 



2.3 Additional Properties 

2.3.1 UandM* 

We derive the host galaxy luminosity, L*, by subtracting the es- 
timated net AGN luminosity from the total observed luminosity. 
For the total observed luminosity we use the SDSS cMODEL flujtj 
(Abazajian et al. 2004) in the z-band, which is a linear sum of a 
de Vaucouleurs model and an exponential model fit to the image, 
and is the best suited model to account for both the galaxy and the 
nuclear light. The z-band is chosen since it is the reddest SDSS 
band, therefore it has the highest host to AGN contrast. It also has 
the smallest dispersion in the ratio of host mass to host light. We 
estimate the net AGN luminosity at the z-band, L^gn; z-band> from 
LbHa (Paper I). 

To convert to A/*, we use the the median mass to light ra- 
tio of (M/L) = 10°- 25 (M/L) Q , found by Kauffmann et al. (2003b, 
Fig. 14 there) for galaxies with L* in the range spanned by the Tl 
sample. This method has an error of ~ 50%. Since the mass to light 
ratio is similar for the SDSS-/ and z-bandfl and all Tl objects have 



1 Not available for seven objects. They are disregarded when M, is used. 

2 The median (M/L) is larger in the /-band than in the z-band by ~ 0.05 
dex, but the normalization is lower by 0.03 dex (Kauffmann et al. 2003b). 



z < 0.3, we do not need to apply a K-correction to L» when calcu- 
lating Mif. The individual M* of the Tl sample objects are listed in 
Table 1. 

An additional source of error is the scatter in the ratio of 
^AGN; z-band to A)Ha- We assume this scatter equals the scatter in 
the relation between L total . 5100A an d A)H/3 °f 0.5 < z < 0.7 SDSS 
quasars - the lower z limit ensures the quasars are luminous and 
host contribution to the continuum is minimal, while the upper z 
limit ensures H/3 fully appears in the spectrum. Using the L(,H/S 
and £< tota j. 5100A vaiues from Shen et al. (201 1), we find a scatter of 
0.2 dex. This scatter implies that in the 9% of the Tl objects with 
implied Lagn /L* > 1 , the true may be overestimated by a factor 
of more than 2, while in the 3% of Tls with implied Lagn/^-* > 3, 
the true L* may be underestimated by a factor of more than 2. In 
0.5% of the objects, the implied L„ is negative. We address these 
limitations when utilizing L* and below. 

2.3.2 L[jy and a ox 

We derive the L uv (= vL v (1528A)) and L x (= vL v (2 keV)) of 
the Tl AGN, from the GALEX (Martin et al. 2005) and ROSAT 
(Voges et al. 1999) surveys. GALEX observed 89% of the Tls, 
and detected 93% of them. ROSAT observed the entire sky, and de- 
tected 43% of the Tls. The derivation of the luminosities is detailed 
in Paper II. Table 1 lists L uv and Oo X = —0.42 x logLuv/Lx, the 
slope of the interpolated power law between the UV and the X-ray. 



2.4 The T2 sample 

We compare our results to the Brinchmann et al. (2004) type 2 AGN 
sample, which was derived from the SDSS galaxy survey, using the 
emission lines measurement of the MPA/JHU grourjj The type 2 
AGN were selected by requiring S/N > 3 in the [O III], H/3, [N II], 
and Ha narrow emission lines, and being above the KeOl 'extreme 
starburst' line in the BPT-[N II] panel. We use all type 2 objects that 
appear in our parent sample (following the S/N > 10 and bad pixel 
cuts, §2.1), excluding the 454 objects which enter the Tl sample, 
as they show broad Ha emission. We name these 13 705 objects as 
the T2 sample. 

The MPA/JHU group modeled the stellar absorption features 
using the Bruzual and Chariot (2003) stellar library. We use a sim- 
pler technique in the Tl sample, based on the Yip et al. (2004) ESa, 
due to possible degeneracies of different stellar components with 
the unobscured AGN continuum (see §2.2 in Paper I). In order to 
understand the effect of the different stellar modeling techniques on 

3 Available at http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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the measured narrow line ratios, and the effect of other differences 
in the fitting procedure, we run our fitting algorithm on 700 spec- 
tra from the T2 samplfl Then, we compare the narrow line ratios 
we measure on these T2s with those published by the MPA/JHU 
group. 

In these 700 type 2s, our algorithm gives [O IIl]/H/3 ratios 
which are on average 0.08 dex larger than the ratios measured by 
MPA/JHU, with a dispersion of 0. 1 1 dex. Our [N Il]/Ha, [S Il]/H0f 
and [O I]/Ha measurements are on average 0.04, 0. 1 and 0.07 dex 
larger than MPA/JHU, with dispersions of 0.09, 0.1, and 0.1 dex. 
The offsets in the narrow line ratios are mainly due to offsets in the 
measured flux of the narrow Ha and Hj8 lines (mean offset —0.09 
dex each), which could imply that we did not fully correct for the 
stellar Balmer absorption features. Therefore, to minimize offsets 
originating from measurement issues, we hitherto decrease the nar- 
row line ratios we measure in the Tl sample objects by these mean 
offsets. 



3 THE BALDWIN EFFECT OF THE NARROW LINES 

In Figure 1, we present the ratio of the narrow lines luminosity 
L nHj3> %n], £<[Sii]> and L [0l] , with L bHa as a function of L bHa . 
Black dots mark objects with robust measurements, while gray 
markers indicate the less robust values. A trend of decreasing NLR 
to BLR luminosity ratio with increasing L\^ a is clearly seen for all 
lines. Similar relations were found in Paper II for [O III] and Ha. 
The least-squares best fit power law slope and associated dispersion 
are noted in the panels. The relations found in Paper II for [O III] 



and Ha are L nH « ~ i^ H 6 > = 0.37 and L [Q m] <* L^ita' a ; 



: 0.36. 



The formal error on all slopes is ~ 0.01. We treat L^ea, which is 
used to select the Tl sample, as the independent variable. 

Since LbHa x A>ol (Paper I), the observed trends represent a 
Baldwin effect (Baldwin 1977) for the narrow lines. Note that the 
different slopes found above imply some trends in the mean posi- 
tions with luminosity of the Tl objects in the BPT plots, as shown 
below. 
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3.1 Less robust values 

In all four panels of Fig. 1, most upper limits fall within the distri- 
bution of the general population. The objects in which the narrow 
lines are fit with an [O IIl]-like profile (§2.2.3) are located at the 
high-LbHa end of the sample. As noted in §2.2.3, the deblending of 
the Balmer lines and [N II] from the broad lines may be inaccurate 
in these objects. Indeed, the L[n n] /A>Ha values of these objects are 
offset to lower values then the general trend. When excluding these 
objects, we find L nH| g oc L b '^ and L[n n] x ^Hcr '- e - a s™ 1 ^ H/3 
slope and a [N II] slope higher by 0.06 compared to when using all 
objects. 



3.2 Comparison with previous studies 

Croom et al. (2002) compared the narrow [O III], [O II], [Ne III], 
and [Ne V] line luminosities with the absolute B magnitudes of 2dF 
and 6dF quasars (Croom et al. 2001). For a direct comparison with 
our results we subtract the slope they found for each line with the 
slope they found for L(,H/i • Comparing the narrow lines to the broad 



The first 700 objects, sorted by right ascension. 



Figure 1. The distribution of L nalrow line /LbHa versus LbHa for narrow lines 
analyzed in the BPT plots below. Robust narrow line measurements are 
marked by black dots. Profiles of [N II] and H/3 based on [O III] (§2.2.3) are 
marked by gray squares. Upper limits (non detections, §2.2.2), and lower 
limits (stellar absorption for H/3 only, §2.2.4) are marked by the appropri- 
ate arrows. The slope of the best fitting power laws (black lines) and the 
associated dispersion are noted. For comparison, the relations found in Pa- 
per II, for [O III] (index = 0.72) and for Ha (index = 0.66), are shown as 
gray lines. The tendency of increasing NLR / BLR ratio with decreasing 
AiHa, found in Paper II, is seen in all lines. Since LbHa 05 L(, i (Paper I), 
the observed trends represent the Baldwin effect of the narrow lines. When 
excluding objects marked by squares, the [N II] Baldwin slope increases to 
0.60. Note that [S II] and [N II] have steeper slopes than Ha, which imply 
a shift in the mean positions in the BPT plots with luminosity. 



H/3 also avoids the bias created by host contamination of the con- 
tinuum (Paper II). Their implied slopes are 0.86, 0.49, 0.58 and 
0.74 for [O III], [O II], [Ne V], and [Ne III], respectively. All their 
narrow lines show a Baldwin effect, as found here. Their [O III] 
slope of 0.86 is steeper then our slope of 0.72, while their [O II] 
slope of 0.49 is flatter than any of the slopes measured here. 
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Mid-IR narrow lines also show Baldwin effects (Honig et al. 
2008, Keremedjiev et al. 2009). 



4 THE BPT POSITIONS OF THE Tl AGN 

Figure 2 presents the BPT positions of the 3 175 Tl AGN, plotted 
over the SDSS narrow line galaxies (Fig. 1 from Ke06). Classifica- 
tion lines are from KeOl, Ka03, Ke06 and Ho et al. (1997a, here- 
after Ho97). The classification of each Tl object in each panel is 
listed in Table 2. 

The T2 AGN reside, by definition, above the KeOl line in the 
BPT-[N II] panel. However, only 80% of the Tl objects reside in 
the AGN regime, 15% are classified as composite and 5% as SF. 
We stress again that all Tl AGN are clearly powered by accre- 
tion onto a massive black hole, as indicated by the detection of a 
broad Ha. Thus, the SDSS type 2 AGN sample is likely only 80% 
complete. Including composites will increases the completeness to 
95%, but may include a significant number of objects which are not 
true AGN. 

We note that the narrow line measurements of two-thirds of 
the Tls which reside in the SF region are poorly constrained. Thus, 
with higher quality spectra the true AGN fraction with SF narrow 
line ratios may therefore be as low as 2%. In comparison, only 18% 
of the Tls classified as composites and 17% of the Tls classified 
as AGN' have poorly-constrained narrow line measurements. 

The fraction of Tl which reside outside the AGN region in 
the BPT-[0 I] panel is 18%, and in the BPT-[S II] panel it reaches 
29%. The SDSS spectra are taken with a 3" fiber, which can include 
a significant fraction of the host galaxy emission. Below we study 
some indications that the offset from the AGN region in the BPT 
plots indeed results from host contamination. 

Fig. 1 also shows that a sizable fraction of the Tl sample oc- 
cupies a new region in the BPT panels, with [O III]/H]3= 5 — 10, 
[N Il]/Ha= 0. 1-0.3, and [S Il]/Ha = 0. 1-0.3. These objects have 
no counterpart in the narrow line sample. Specifically, 10% of Tls 
with [O III] /H/3 > 5 have [N Il]/Ha < 0.3, compared to only 0.8% 
of the T2 sample. Below, we study the range of AGN and host 
properties at which these line ratios are dominant, and discuss their 
physical origin. 

In the [O I] panel, 190 objects in the Tl sample are classified 
as LINER Is. The vast majority (179) of them appear to the right 
of the Ho97 line, where the 'bona-fide' LINERs reside. This result 
is consistent with the strong drop in broad Ha detection across this 
line (Ho et al. 1997b, Ho 2008, Wang et al. 2009). 

4.1 BPT positions of Tl AGN, by LbHa* L/LEdd and Lagn /L* 

In this section we utilize the large size of the Tl sample, and ex- 
plore their positions within the BPT plots when the sample is cut 
based on various AGN and host properties. We identify some qual- 
itative trends, which are further explored in the following sections. 

Figure 3 presents the BPT positions of the individual Tl 
objects as a function of LbHa> which is a measure of L^i (= 
130 x LbHa, Paper I). As in Fig. 2, the positions of the Tls 
are plotted over the SDSS narrow line galaxies from Ke06. At 
logL bHcf = 40.7 (L bo i = 42.8), Tl AGN largely overlap the nar- 
row line sample. With increasing luminosity, the Tl AGN shift to 
lower [N Il]/Ha, lower [S Il]/Ha values, slightly lower [O I]/Ha, 
and higher [O IIl]/Hj8 values, as expected from the different lumi- 
nosity trends of the different lines (Fig. 1). At quasar luminosities 
(log LbHa > 43 or logLb i > 45), the Tl distribution is distinct 



from the type 2 distribution in the [S Il]-panel, and is offset in the 
[N Il]-panel to lower values. Also, the fraction of AGN which re- 
side below the KeOl line decreases with increasing LbHa- 

The fraction of poorly-constrained objects increases with 
LbHa (Fig- 1)^ due to the decrease in the relative strengths of the 
narrow lines (Paper II and Fig. 1). Therefore, one may wonder 
whether this trend with luminosity is not simply due to the limi- 
tations of the deblending algorithm. In appendix A we verify the 
observed trend using high quality mean spectra. 

Figure 4 shows the BPT positions of the Tl objects, now sub- 
divided by L/L Edd (derived from LbHa and Av, using eq. 3 in Paper 
I). Several trends in the BPT position are apparent. With increasing 
L/L Edd , an increasing fraction of Tls have highs [O III]/H]3, low 
[N Il]/Ha, and low [S Il]/Ha, as found with increasing LbHa m 
Fig. 3. The decrease in [O I]/Ha with L/L Edd is more pronounced 
than in Fig. 3: the [O I]/Ha span mainly 0.1 - 0.3 at low L/L Edd , 
compared to 0.03 - 0.1 at high L/L Edd . The fraction of LINERs 
in the BPT-[0 I] panel strongly decreases with increasing L/L Edd , 
from 29% at log L/L Edd = -2.5, to 6% at log L/L Edd = -1.8 and 
3% at log L/L Edd = -1.2 and -0.6. 

Figure 5 is similar to Figs. 3 and 4, with different rows des- 
ignating different bins in Lagn A* , measured at the SDSS-z band 
(see §2.3.1). We note the division of objects between the two high 
Lagn/L* bins is not robust in objects with Lagn/^* > 1> due to 
the possible error in the determination of L*. The 24 objects with a 
negative implied L* appear in the Lagn /£* ^ 2 bin. With increas- 
ing Lagn/^*, [Nll]/Ha and [S Il]/Ha decrease, as found with 
increasing LbHa m Fig- 3, and with increasing L/L Edd in Fig. 4. 
The Tl sample overlaps the type 2 sample in host dominated ob- 
jects, and is distinct from the type 2 distribution in AGN domi- 
nated objects. Also, the composite fraction decreases from 22% at 
Laqn/U = 0.04 to 6% at L ACN /U > 2. 

To satisfy the curious reader, the mass dependencies are ex- 
plored in Appendix B, where we plot the BPT positions of the Tl 
objects, subdivided by Mbh and M*. 

4.2 Comparison with previous type 1 samples 

Some earlier studies have measured the BPT positions of smaller 
samples of type 1 AGN. Dietrich et al. (2005), compiled the BPT 
positions of earlier Seyfert 1 galaxies in the literature (Cohen 1983, 
Ho97), and found that Seyfert Is largely overlap Seyfert 2s (Fig. 
7 there). This result is consistent with our finding (Fig. 3, upper 
panel), that type 1 AGN at low luminosity overlap the positions of 
type 2 AGN. 

Greene & Ho (2007) inspected the narrow line ratio of 229 
SDSS type 1 AGN, selected based on the detection of a broad Ha, 
similar to Tl, but required to have Mbh < 2 x 10 6 M Q . In their 
sample, 39% of the objects are classified as Composites or SFs in 
the BPT-[N II] panel, versus only 20% (Fig. 2) in our sample. How- 
ever, when we restrict the Tl sample to M B h < 2 x 10 6 M© (Fig. 
Bl, upper panel), the fraction increases to 36%, consistent with the 
Greene & Ho (2007) result. In a followup paper (Xiao et al. 2012), 
they compared the narrow line ratios based on the SDSS spectra 
with ratios based on spectra from a smaller aperture. The fraction 
of Composites / SFs decreased to 18%, indicating that extended 
emission from SF in the host galaxy shifts the BPT position into 
the composite region at low M B h- Below, we provide further evi- 
dence that this effect applies also to composites at higher Mbh- 

Another prominent feature in the Greene & Ho low Mbh 
sample are objects with high [O III]/H|3, low [S Il]/Ha and low 
[N Il]/Ha. Of their type l's with [O IIl]/Hj3 > 5, 10% have 
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Figure 2. The BPT positions of the 3 175 Tl AGN (solid contours), compared to the narrow line SDSS galaxies (background gray pixels, Fig. 1 in Ke06). Red 
solid lines are the KeOl extreme starburst lines, while the dashed lines mark the separation between 'SF' and 'Composites' in the [N II] panel (Ka()3), and 
between 'Seyferts' and 'LINERs' in the [S II] and [O I] panels (Ke06). The dash-dotted line in the [O I] panel further divides the LINER group into bona-fide 
LINERs and 'transition' objects (Ho97). The contours encircle regions with 10, 25, 60 and 150 Tl objects per 0.2x0.2 dex 2 bin, respectively. Beyond the outer 
contour, Tl AGN are marked as in Fig. 1, with upper / lower limits on either of the emission lines denoted by an appropriate arrow. In the [N II] panel, 20% 
of the Tl objects are below the KeOl line, and would not be defined as AGN. Of these, 15% are defined as Composites, and 5% as SF galaxies. In the [O I] 
panel, 190 objects in the sample are LINER Is, of which 94% (179/190) appear to the right of the Ho97 line. Note that the Tl extends to higher [O III]/Hj6, 
and lower [N II]/Ha and [S II]/Ha, compared to the narrow line galaxies and AGN. 



[N Il]/Ha < 0.3, as found for the entire Tl sample above. These 
narrow line ratios are also observed in low Mbh type 2 samples 
(Barth et al. 2008), but are clearly missing from type 2s at higher 
luminosity (Fig. 2). It therefore seems that type 1 and type 2s have 
similar ratios at low Mbh, but become distinct at higher Mbh- We 
examine the reason for this difference below. 

Winter et al. (2010) published the BPT positions of a hard X- 
ray selected AGN sample, of which they identified 33 objects as 
broad line AGN0 Their mean log L[o m) ls a factor of three higher 
than the mean in the Tl sample. Five of their type ls have [S Il]/Ho! 
< 0.1, and two have [N Il]/Ho! < 0.1, which are not seen in their 
type 2 sample. These ratios are seen in the Tl sample, but not in 
the SDSS narrow line sample. 

Buttiglione et al. (2010) measured the BPT positions of a ra- 
dio selected AGN sample. They show a clear trend of decreasing 
[N Il]/Ho! and [S Il]/Ha, and increasing [O IIl]/H/3, with increas- 
ing [O III] luminosity (Fig. 1 there). Their trend is equivalent to 
the trend seen in Fig. 3 with L\,u a . A similar trend can be seen in 
Fig. 4 of Wang & Wei (2010), who measured the BPT ratios of 
Seyferts 1.8s and 1.9s with ROSAT detections. They found that ob- 
jects in which the AGN dominates the continuum are offset to lower 
[N II]/Hgi than objects in which the continuum is host dominated, 
as can be seen here in Fig. 5. 



5 They excluded Syl.8s and Syl.9s, which selects against low luminosity 
type 1 AGN (Paper II). 



5 THE NATURE OF Tl AGN WITH LOW [N n]/Hce AND 
LOW [S n]/Ha 

Figs. 2-5 show that the distributions of the narrow line ratios of the 
Tl objects extend to values which are not seen in the SDSS nar- 
row line sample, in particular at high L\, \, high LjL^&A, and high 
^AGN / L* . These non overlapping objects have [O IIl]/H/3 similar 
to type 2s, but lower [N Il]/Ha and [S Il]/Ha. In Appendix A we 
show the offset ratios are not an NLR/BLR deblending artifact. 
Why are these ratios absent from type 2 samples? what is the phys- 
ical mechanism that produces them? 

5.1 The difference between the Tl and T2 samples 

A difference between type 1 and type 2 AGN can be either a failure 
of the unified model, an orientation-related effect, or simply due to 
different selection criteria used for creating the two samples. Here, 
we compare the Tl and T2 (§2.4) samples, and show that selec- 
tion effects are behind the differences observed in Figs. 2-5. To 
avoid significant NLR contamination by star formation in the host, 
which decreases the narrow line ratios to the Composite and SF re- 
gions of the BPT plots, we only use the 1 691 Tls and 4042 T2s 
with [O IIl]/H/3 > 5. This selection criterion is independent of the 
offset quantities - [N Il]/Ho: and [S Il]/Ha. A comparison of the 
Composites and SFs is made in §7 below. 

A major difference in the Tl versus T2 selection criteria, is 
that the Tl sample includes also point sources, and is not selected 
purely from extended objects. Thus, the Tl sample can extend to 
Lagn/L* values larger than possible in the T2 sample. Another 
related systematic difference is the distribution of M* values, as the 
T2 objects are selected by L„, while in the Tl point sources L* can 
be arbitrarily small. The distribution of L„ values is interesting as 
M* was found to correlate with the [N Il]/Ha ratio, via the M* — Z 
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Figure 3. The dependence of the BPT position on the AGN luminosity. The Tl sample markers are as in Fig. 1. The narrow line background and dividing 
lines are described in Fig. 2. Each row presents Tl AGN from a given decade-wide bin in Z-bHa- The mean LbHa in each bin are noted (in erg s ), At 
logZ.t,Ha = 40.7 (Lb„i = 42.8), Tl AGN overlap the narrow line sample. With increasing luminosity, the mean [N II]/Hgi and [S II]/Ha decrease. At quasar 
luminosities (log LbHa ^ 43 or logLbol > 45), a large fraction of the Tl AGN occupy a region in the [N II] and [S II] panels which is distinct from the type 2 
distribution. 



relation of galaxies, and the dependence of [N Il]/Ha on Z (Groves 
et al. 2006, see below). 

In Figure 6, we therefore plot contours of the distribution of 
T 1 s and T2s in the vs. L^gn / ^* plane. The T 1 sample is divided 
according to the two SDSS surveys from which it is derived, those 
selected from the SDSS galaxy survey, and the point-sources from 
the SDSS quasar survey. We note that at Lagn/-^* > 1, £<* (and 
M t ) can be significantly overestimated (§2.3.1), therefore the true 



^-AGN /L* may be higher and the true M» may be lower than plotted. 
In the 17% of Tl point sources with L^Gn/^* > 3, L t is not well- 
constrained either from below or from above, therefore we abruptly 
cut the distribution. This cut does not affect the conclusions below. 
In the T2 sample, Z-agn ls derived from L[o m] (Paper II), and L„ 
and M» are derived from the observed SDSS z-band luminosity, 
using the same M/L as in Tls. 

The T2s are all selected from the SDSS galaxy survey. The 



Type 1 AGN at low z. III. The optical narrow line ratios 9 




log [Nil] /Ha log [Sn]/Hce log [Orj/Hce 



Figure 4. As Fig. 3, for the dependence of the BPT position on L/Lem- Each row presents Tl AGN for a given L/L^dd bin (0.75 dex wide). The mean L/LecM 
in each bin is noted in the lower left corners. With increasing L/Lfcdd, 'he Tl AGN move to the left in all BPT panels. Also, at log L/^Edd = —2.5, 29% of the 
Tl sample are in the LINER region of the [O I] panel, compared to 6% in the entire Tl sample. 



distribution of the Tls from the same survey overlaps well the dis- 
tribution of T2s. The Tl point sources however are clearly off- 
set to higher Lagn/^*! and constitute an AGN population which 
does not appear in the T2 sample. As Fig. 5 shows, the Tls be- 
come offset from T2s at the higher Lagn/^* values. Thus, the 
apparent differences between the Tl and T2 BPT positions re- 
flects their different Lagn/^* values, which controls the BPT po- 
sitions. Figure 7 explores this effect more quantitatively. The Tls 
and T2s are binned by the [N Il]/Ha values, and the distributions 
of the different bins in the vs. Lagn/^* plane is presented. 



In the Tl sample, M„ decreases and Lagn/^* increases with de- 
creasing [N Il]/Ha. The T2 sample shows a similar trend, but it 
does not extend to the high Lagn/^* occupied by objects with 
[Nll]/Ha < 0.2. In fact, already T2s with [N Il]/Ha < 0.6 are 
rare, constituting only 8% of the T2 sample, compared to the 39% 
of Tls that have [N Il]/Ha < 0.6. Therefore, the T2 objects do not 
extend to the low [N Il]/Ho: values, seen in the Tl sample, as these 
values occur at high Lagn/^* values, which the T2 objects cannot 
have by their selection. 

A similar analysis using the [S Il]/Ho! values, instead of 




Figure 5. As Fig. 3, for the dependence of the BPT position on Lagn/£*- Each row presents Tl AGN from a given Lagn/^** bin (0.7 dex wide). The Lagn/^* 
is calculated at the SDSS z-band, and the mean values are noted in the lower left corners. The Tl AGN move to lower [N II]/Hce and [S II]/Ha with increasing 
Lagn/L,, similar to the trend with L^n a seen in Fig. 3. The composite fraction decreases from 22% at L^qn/L, = 0.03 to 6% at Lagn/£* ^ 2. 



[N Il]/Ha, demonstrates that the low [S Il]/Ha values seen in the 
Tl sample at high L^gn/^* (Fig- 5) are absent from the T2 sample 
for the same reason. 



5.2 The quantitative change in [N n]/Ha and [S n]/Ha 

After addressing why the low [N Il]/Ha and [S Il]/Ho! ratios seen 
in the Tl sample are absent from the T2 sample, we now quantify 



the dependence of [N Il]/Ha and [S Il]/Ha on the observed AGN 
and host properties. 

As can be seen in Figs. 3, 5 and 7, [N Il]/Ha decreases 
both with increasing Lagn an d with decreasing L*. Below we ex- 
plore the two effects independently. We bin the Tl objects with 
[O IIl]/H/3 > 5 (to avoid host contamination) based on M t and 
L), i (= 130 x LbHa) in me following manner. The objects are 
sorted by ibol an d divided into four equal size groups. Each of these 
groups is then sorted by M„, and again divided into four equal size 
groups. This ensures similar statistical errors in all bins. We disre- 
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Figure 6. The distribution of the Tl and T2 AGN, with [0 IIl]/H/3 > 5, in 
the M t versus L^gn/L* plane. The distributions are denoted by two con- 
tour lines, which indicate the number of objects per 0.3 x 0.3 dec 2 , 1% of 
the size of the group for the outer contour, and 2.5% for the inner contour. 
The three groups marked are Tl selected as point sources (thin solid), Tl 
selected to have extended morphology (thick solid), and T2s which are all 
selected to have an extended morphology (dotted). At Lagn/L* > 1, L.,_ 
(and M„) may be overestimated (§2.3.1). In the 17% of Tl point sources 
with Lagn/L. > 3, L t is not well-constrained either from below or from 
above, therefore we abruptly cut the distribution. The extended Tl distribu- 
tion overlaps the T2 distribution, as expected from their common selection 
criteria (and AGN unification). The Tl point sources, selected by their non- 
stellar colors, are offset to higher Lagn/L, than AGN with an extended 
morphology, and constitute an AGN population which does not appear in 
the T2 sample. 



gard the 12% of the objects with L^gn/L* > 1, in which the M„ 
measured has a large error (§2.3.1). 

Figure 8 presents the derived relations of the mean values of 
[N Il]/Ha, [S Il]/Ha, [N II]/[S II], and L nlia /L blia as a function of 
M«, for different Lbol- Error bars denote the error in the mean. The 
upper panel shows that the mean [N Il]/Ho! increases with M r at a 
fixed L(, i, an increase of ~ 0.3 dex over ~ 0.7 dex in M*. Also, the 
mean [N Il]/Ho! decreases with L D ol> at a fixed M*, a decrease of 
~ 0.3 dex over ~ 1.7 decades in L\, \. The second panel shows that 
[S Il]/Ho! increases only slightly with M- t at a fixed L\, \. However, 
the mean [S Il]/Ho! decreases with Lb i at a fixed M„, similar to 
the decrease in [N Il]/Ha with L DO i in the upper panel. The rela- 



Figure 7. The distribution of Tl and T2 AGN samples with [O III]/H/3 > 
5 in the M, versus Lagn/L» plane. Each sample is grouped into bins in 
[N II]/Ha. The contour line of each bin is equivalent to the inner contour 
of each group in Fig. 6. The [N II]/Ha range of each bin is noted near 
the contour, and coded by color. Tl contours are solid, and T2 contours 
are dotted. In the Tl sample, the distribution moves to lower M* and higher 
Lagn/L, with decreasing [N II]/Hce. T2s follow a similar trend, but they do 
not extend to the high Lagn/L» occupied by objects with [N II]/Ha < 0.2 
(Fig. 6). Thus, the absence of [N II]/Ha < 0.2 AGN from the T2 sample 
just reflects its selection against high Lagn/L* objects, which are observed 
to have a low [N ll]/Ha < 0.2. 

tive trends of [N II] and [S II] are most apparent in the third panel. 
Clearly, [N II]/[S II] strongly increases with M*, and is almost in- 
dependent of L DO i, w i m only a small decrease at the highest Lb i 
bin, at low M* . 

What are the physical mechanisms behind these trends in 
[N Il]/Ha and [S Il]/Ha? The similarity of the behavior of 
[N Il]/Ha and [S Il]/Ha vs. L DO i, in contrast to the different be- 
havior vs. Ms, suggests there are two distinct mechanisms at play. 
We address them separately below. 

5.3 The trend with M* 

5.3.1 M„ vs. Z NLR 

An increase of [N Il]/Ha with M, has been observed in type 2 AGN 
by Groves et al. (2006), qualitatively similar to the trend we see 
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Figure 8. The dependence of the mean line ratios of the Tl objects (with 
[O III]/Hj8 > 5) on Lb i and M t . Given Lbol bins are connected by solid lines 
(mean log i^ol noted). Error bars denote the uncertainty in the mean value. 
Each of the 16 bins in the top and bottom panels has 92 objects, while each 
bin in the middle panels has 76 objects. Top panel The mean [N II]/Ha 
increases with M t at a fixed Lbol, and decreases with increasing Lbol at a 
fixedM,. Second panel The mean [S II]/Ha decreases with increasing Lbol 
at a fixed AL,. At a fixed Lboi, [S II]/Ha only slightly increases with M,. 
Third panel The mean [N II]/[S II] is determined by M„ and is almost 
independent of Lbol ■ The three trends with M„ shown above likely reflect an 
increase of Zjmlr with M». Bottom panel The mean L n Ha/ibHa decreases 
with Lboi an d is largely independent of M t . The LnHa/A>Ha is 3 measure 
of the NLR covering factor (Paper II). Therefore, the decrease in [N II]/Ha 
and [S II]/Ha with Lboi, at a fix M», is likely related to the decrease in NLR 
covering factor with increasing Lbol ■ 



in Tl AGN (Fig. 8). As mentioned above, Groves et al. claimed 
this trend originates from the M* — Z relation found in quiescent 
galaxies. 

A dependence of [N Il]/Ha on Znlr is expected since Nitro- 
gen is a secondary nucleosynthesis product, and hence its abun- 
dance increases as Z 2 for Z > 0.5 Z Q (e.g. van Zee et al. 1998). 
This secondary nature also implies that the relative abundances of 
N and S in the NLR increases °c Znlr- Fig. 8 shows clearly that the 
[N II]/[S II] increases with M„, consistent with a change in Znlr 
with M* . 

In the photoionization models of Groves et al. (2004, 
Fig. 19 there), [N Il]/Ha increases roughly like Z^ R , com- 
pared to [S Il]/Ha which increases roughly like Z^r, supporting 
[N II]/[S II] cx (N/S) NLR °c Z NLR . Hence, the [N II]/[S II] panel of 
Fig. 8 implies a change of 0.2 dex in Znlr over 0.8 dex in M*. 
For comparison, Groves et al. (2006) find a change of 0.2-0.3 dex 
in Znlr over two decades in M,. We speculate that our steeper 
— Znlr relation is because that by using [N II]/[S II], we fac- 
tored out the decrease of [N Il]/Ho: with L bol (upper panel in Fig. 
8), which flattens the relation found by Groves et al. 

In inactive galaxies, Tremonti et al. (2004) found an increase 
of ~ 0. 1 dex in the mean gas phase Z at 10.4 < log M* < 11.2, half 
the increase of 0.2 dex in Znlr estimated here. 



5.3.2 Znlr vs. Zblr 

Are there any additional differences in the spectra of objects with 
high and low [N Il]/Ha? Figure 9 compares the mean spectra of 
objects with [N Il]/Ha < 0.2 and objects with [N Il]/Ho: > 0.6. 
To avoid other known trends, and isolate only [N Il]/Ha related 
trends, we match each of the Tls with [N Il]/Ha < 0.2 with a Tl 
that has [N Il]/Ha > 0.6 with the same L\, \ up to 0.1 dex, and the 
same Av up to 0.05 dex. Matching by L(, ] ensures we are freezing 
the L(, i-related effect seen in Fig. 8, while matching also by Av 
indicates we are freezing also Mbh (via eq. 2 in Paper I) and L/LEdd 
and the host of spectral properties related to it (e.g. Boroson & 
Green 1992). Of the 98 Tls with [N Il]/Ha < 0.2, 79 have such 
matches. 

The mean spectra of the two groups of objects are calculated 
by geometrically averaging luminosity densities of spectrum pixels 
with the same restframe wavelength A, rounded to 10~ 4 in logA. 
The bottom spectrum is the difference between the two composite 
spectra, and the insets zoom in on the areas delimited by the dashed 
lines. The most striking feature of the residual is the strong BLR 
Fe II multiplets at ~ 4600A and ~ 5300A. 

The luminosity of the optical Fe II multiplets is expected to 
increase with iron column density, and therefore with Zblr, to a 
power of 0.8-0.9 (Vemer et al. 2003, Baldwin et al. 2004, Shields 
et al. 2010). Thus, Fig. 9 provides interesting evidence that the BLR 
metallicity, Zblr, is related to Znlr- There is a well-known relation 
between the Fe II equivalent width and L/LEdd (Boroson & Green 
1992), but since the two composites are matched in L/L^dd, this 
effect should not be present. 

Shields et al. (2010) found that when binning by 
L(Fe Il)/Lt,H/j, [N II]/[S II] increases by a factor of two for 
an increase of a factor of ten in L(Fe II) /L bli p . They concluded 
that the Fe II strength increases with Znlr, but the dispersion in 
Fe II is not dominated by Znlr- In Fig. 9, the composite spectra 
differ by a factor of ~ 3 in [N II]/[S II], and a factor ~ 2 in L(Fe II), 
suggesting that for a constant L/LEdd, Znlr an d Zblr change in 
unison. 
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Figure 9. The difference between mean spectra of Tl AGN with high (> 0.6) and low (< 0.2) [N II]/Ha values. The two samples are selected to match in 
LbHa an d Av, and thus in L/Lgdd as well, Insets zoom in on the dashed rectangles in the difference spectrum. The high [N II]/Ha spectrum also has stronger 
BLR Fe II multiplets at ~ 4600A and ~ 5300A, which are expected to increase with Zblr- Therefore, at a fixed L\, a \ and L/Lem, Znlr and Zblr appear to 
be related. The high [N II]/Hce spectrum also shows a larger M t , evident from the red slope and the stellar absorption features at 3934A and 6500A in the 
residual. Thus, Zblr also appears to be at least partly driven by the host M t . 



The mean log M t of the low and high [N Il]/Ha composite 
spectra are 10.5 and 10.9, respectively. This difference in can 
be seen in the residual spectrum, which has a red optical slope, a 
[Ca II] K A3934 absorption feature, and a stellar absorption blend 
at 6500A. The two groups are selected to have the same mean Mbh, 
and should thus have similar mean bulge mass (Magorrian et al 
1998). The different measured mean M* values of the two groups 
should therefore reflect differences in the mean disk masses, where 
the higher metallicity group has a higher disk/bulge mass ratio. 

Hamann & Ferland (1993, 1999) found that in quasars, Zblr 
(derived from the NV / CIV ratio) increases with Z^ol- They spec- 
ulated that the increase in Zblr with ibol 1S probably due to the 
increase of Zblr with M.* > an d the strong relation between M* and 
L(, ] in the quasar samples they used. Here, we confirm their claim 
by showing Zblr increases with M* directly. At a given M*, the 
mean [N II]/[S II] remains constant with L(, i (Fig- 8), and only 
slightly decreases at the highest L\, \. Furthermore, this slight de- 
crease could be due to improper deblending of the [N II] profile 
from the broad Ha profile at high L^ea ( see Fig- 1)- Therefore, we 
find no evidence for a direct Z — L^oi trend. 



5.4 The trend with Lt, i 

As mentioned above, the decrease of [N Il]/Ha with L\, Q \ is proba- 
bly not a Znlr effect. Therefore, deriving Znlr in quasars from 



[N Il]/Ha, using the calibration of lower luminosity AGN (e.g. 
Husemann et al. 201 1), will underestimate Znlr- 

A hint on the physical source of the trend with L\, a \ can be 
seen in the bottom panel of Fig. 8. The i n Ha/^bHa decreases with 
Lbol> as found in Paper II, but is only weakly dependent of Af», if at 
all. Therefore, it seems that the decrease in [N Il]/Ha and [S Il]/Ha 
with Lb i, at a given M- t , is related to the decrease in L n Ha /A>Ha ■ In 
Paper II, we presented evidence that the L n ^ a /L\^ a vs. i^ol trend 
is due to a decrease in the NLR covering factor with L^ol- The phys- 
ical connection between the decrease in [S Il]/Ha and [N Il]/Ha 
and the decrease in NLR covering factor remains an open question. 



6 LINERS 

Ke06 found that at a fixed L/L^^, the difference between host 
properties of Seyferts and LINERs disappear. Their conclusion was 
that the observed difference in host properties between Seyferts and 
LINERs is only a secondary effect, which results from their differ- 
ence in L/Lgdd (Ho 2002, Ke06). Here, we show that the observed 
large difference between Seyferts and LINERs in terms of the frac- 
tion which shows broad lines (Ho et al. 1997b, Ho 2008), is also a 
secondary effect of their difference in L/L^^, and at a fixed L/L^dd 
the difference disappears. 

Following Ke06, we create subsamples of the Tl and T2 sam- 
ples which include objects classified as AGN in the BPT-[N II] 
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panel, and as either Seyferts or LINERs in the BPT-[S II] and BPT- 
[O I] panels (Fig. 2). We use only objects with consistent BPT-[S II] 
and BPT-[0 I] classifications. We use the bulge stellar dispersion 
o t to derive M Bli in T2s (Giiletkin et al. 2009). We disregard the 
8% of the T2s with surface mass density < 3 x 10 8 Mq kpc~ 2 , in 
which the cx* measured by the SDSS may be overestimated due to 
disk light contamination (Kauffmann et al. 2003c, Heckman et al. 
2004). 

Following the above criteria, the T2 subsample includes 4 938 
Seyfert 2s and 4 292 LINER 2s. The Tl subsample includes 1910 
Seyfert Is and 76 LINER Is. Thus, LINERs constitute 50% of the 
T2 sample, but only 4% of the Tl sample. Our purpose is to further 
understand the origin of this large difference. 

In 44 objects from the LINER 1 group, the classification is 
ambiguous, either due to upper/lower limits on the BPT ratios, or 
because their narrow line ratios are poorly constrained (§2.2.3). We 
address this uncertainty below. The fraction of Seyfert Is with an 
ambiguous classification is negligible. 

In the Tl sample, we derive the L/L E dd from LbHa an d Av, 
using eq. 3 in Paper I. In the T2 sample, we derive L\, \ from L[o m] 
(Paper II), which was derived from the Lp m] - LbHa relation in the 
Tl sample (eq. 3 there). We note in passing that in LINERs, almost 
by definition, L[q hi] /LbHa is expected to be lower (see factor of 
two drop in L[q hi] /^bHa m the lower-left panel of Fig. 6 in Paper 
II), so the true L\, \ in LINERs is likely a bit higher than our derived 

Aral- 

Figure 10 presents the fraction of LINERs 2 out of the T2 sam- 
ple, as a function of L/Lgdd- Seyfert 2s and LINER 2s are cleanly 
separated in L/L E dd, as found by Ke06. At L/L E dd > — 2 all T2s 
are Seyferts, while at logL/L Edd < -4 all T2s are LINERs. Fig. 
1 also presents LINER 1 fractions out of the Tl sample, at differ- 
ent L/LEdd- The uncertainty in the LINER 1 fraction is due to the 
44 Tl objects with an ambiguous LINER classification. At a fixed 
L/L E dd, the fraction of LINER Is is consistent with the fraction 
of LINER 2s. Therefore, the small fraction of LINERs in the Tl 
sample results from the fact that the sample does not extend to low 
enough L/L E dd, where LINERs become the dominant population. 

Why do the Tl and T2 samples differ in their L/Lgdd distribu- 
tion? The fact that all Tls have L/L E dd > 10~ 3 could be a detection 
limit, since low L/L E dd have weak and wide broad Ha features, 
which are hard to distinguish from the stellar continuum (see Fig. 6 
in Paper I). Alternatively, there may be a physical reason for a lack 
of low L/Lgdd type 1 AGN, related to the subject of 'true type 2' 
AGN (see further discussions in Laor 2003; Laor & Davis 201 1). 

The fact that the probability an object is a Seyfert or a LINER 
does not depend on whether the BLR is detected or not, indicates 
the transition from Seyferts to LINER does not affect the BLR. This 
may indicate that the physical difference between these two type of 
objects occurs beyond the BLR, and hence external to the central 
source. I.e., Seyferts and LINERs may differ by the conditions in 
the circumnuclear gas, and not by a different accretion mode. Such 
a scenario implies that the intrinsic UV and X-ray emission of LIN- 
ERs and Seyferts should not be distinct, as found by Maoz et al. 
(2005, 2007). Though, these latter results are disputed (see review 
in Ho 2008). 



7 Tl AGN CLASSIFIED AS COMPOSITES AND SF 

Why do some of the Tl objects display narrow line ratios char- 
acteristic of Composites and SF galaxies? Can such line ratios be 
powered by accretion onto a massive black hole, or does it result 
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Figure 10. The fraction of LINERs in the Tl and T2 samples, as a func- 
tion of L/LEdd- The solid gray line indicates the fraction of LINERs in the 
T2 sample, in 0.25 dec L/Lgdd bins, where L/Lgdd is derived from L[o m] 
and o~*. The fractions of LINERs in the Tl sample are denoted by error 
bars, where L/L Ed d is derived from L bHa and Av. The uncertainty is due 
to Tl objects with an ambiguous classification. At log L/LEdd > -2 all 
T2s are Seyferts, while at logL/L Ed d < -4 all T2s are LINERs. At a fixed 
L/Lndd, the fraction of LINER Is is consistent with the fraction of LINER 
2s within the uncertainties. The low fraction of LINERs in the Tl sample 
(~ 4%) versus the high fraction in the T2 sample (~ 50%) results from the 
difference in the L/Lgdd distribution of the Tl and T2 samples. The lack of 
i/^Edd < 10~ 3 Tls could be due to detection limits, or due to a physical 
absence of low L/L Ed d type 1 AGN. 



from host contamination? The fraction of Tls classified as Com- 
posites increases with L*/Lagn at the SDSS-z band (Fig. 5), which 
suggests a host contamination effect. Below we explore quantita- 
tively the host contamination, based on other indicators, and its re- 
lation to the narrow line ratios. We compare the L nH a/LbHa and 
i-Uv/ibHa °f Composites with those of Tls which fall above the 
KeOl line in the BPT-[N II] panel (hereby called 'pure-AGN'). In 
pure- AGN L nHa and Ltjv correlate with L^a, thus host contribu- 
tion should manifest as higher L n Ha/LbHa and Luv/LbHa due to 
line and continuum emission from the SF regions. 

In Table 4, we list the geometrical mean of LbHa and 
AiHa/A>Ha f° r the Tl AGN classified as pure-AGN, Composites 
and SF. The SF group is divided into 'SF-robust' (32 objects) and 
'SF-non robust' (69 objects), depending on whether their narrow 
line ratios are well-constrained (§2.2.3). This division is to guard 
against systematic uncertainties in the less secure measurements. 
As seen in the lower left panel of Fig. 4, non-robust SFs tend 
to have high L/L E dd, where the NLR is weak and the broad Ha 
is relatively narrow, making the NLR / BLR deblending difficult. 
It is therefore possible that in non-robust SFs broad Balmer flux 
was mistakingly assigned to the narrow Balmer lines, and their SF 
classification is not real. In the Composite and pure-AGN classes 
poorly constrained objects are less abundant (17%), and therefore 
a separate group is not required. 

Since the mean AGN L nHa /LbHa decreases with increasing 
AGN luminosity (Paper II), we compare each classification with a 
pure-AGN matched in LbHa - The matched groups are constructed 
by randomly selecting 1^1 pure-AGN Tl objects with the same 
LbHa ( U P to 0.1 dex), for each Composite or SF (see Table 4). The 
geometrical mean L n Ha/LbHa of the Composites is 0.27, compared 
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BPT-[N II] Classification N L bHa Lai a /L baa L V y/L bHa 



pure- AGN 2303 42.0 0.11 33 

Composites 407 41.8 0.27 54 

L matched pure- AGN 407 41.8 0.12 36 

SF-robust 32 41.6 0.44 ±20% 62 ±20% 

L matched pure- AGN 128 41.6 0.14±10% 40±10% 



have a higher mean Luv /A>Ha> as found earlier, also the highest 
Luv/^-bHa AGN are on average composite in nature. Therefore, 
the excess UV, seen in low luminosity AGN, likely arises from star 
formation in the host, as suggested in Paper I, based on the pure 
AGN SED. 

In addition, AGN within the 'pure- AGN' BPT regime can also 
be affected by host contamination, in particular when getting close 
to the Ke06 line. Narrow emission lines, powered purely by accre- 
tion, likely produces a smaller dispersion than observed in the BPT 
plots. 



8 THE IONIZING SPECTRUM SEEN BY THE NLR 

8.1 <Xox as a measure of the ionizing spectrum slope 

What produces the scatter in the BPT plots? Possible parameters 
are the ionizing spectral slope and the ionization parameter (e.g. 
Groves et al. 2004). Below we test this explanation by exploring 
the dependence of the BPT positions on a ox , the power law slope 
interpolated from Luv and Lx- We use the 752 Tl objects that were 
observed by GALEX and have LbHa > 10 42 ' 5 erg s~ 1 , to avoid host 
contamination of the UV. These Tl objects are divided into bins of 
CHqx with width of 0.25. Figure 12 shows the mean and scatter of 
the BPT-[0 I] positions of the different a ox bins. We use the BPT- 
[O I] panel since it is most sensitive to the ionizing slope (Groves et 
al. 2004). The X-ray detection rates are 77%, 73%, 67% and 29%, 
for the Cfcx = -1.2, -1.4, -1.6 and -1.8 bins, respectively. The UV 
detection rate is 60% for the a ox = — 1.2 bin, and > 95% in the 
other bins. Upper limits are used when a detection is not available, 
so the true a ox of the a ox = —1.8 bin is likely < —1.8, while the 
true a ox of the a ox = -1.2 bin is likely > -1.2. 

For comparison, Fig. 12 also shows the expected BPT-[0 I] 
position for ionizing spectra with different slopes and for different 
ionization parameters, taken from Fig. Id in Groves et al. (2004), 
which assume a density of 1000 cm~ 3 and Z = 2 Z Q . Clearly, the 
observed mean position is independent of the mean observed a ox , 
in sharp contrast with the models which predict a strong depen- 
dence. This discrepancy may either indicate that a ox is not a good 
indicator of the ionization slope at the EUV, or it may indicate that 
the dispersion is produced by another parameter, such as Z. Another 
mechanism that can produce a dispersion in a ox , with no effect on 
the BPT positions, is dust along the line of sight (Paper I). Nor- 
mal dust has a higher absorption at UV wavelengths than at X-ray 
wavelengths. If this dust is located beyond the NLR, the NLR will 
see the intrinsic ionizing spectrum, and the BPT ratios will not be 
affected. 



8.2 The effect of dust on the ionizing spectrum 

Some AGN appear to be dust reddened based on their SED (e.g. 
Richards et al. 2003). In Paper I, we found that the Luv/Arfta dis- 
tribution at the high luminosity end of the Tl sample is at least 
partially due to dust reddening along the line of sight. In particular, 
objects with Luv /LbHa < 30 show a correlation such that redder 
optical slopes go with a decreasing Luv/^bHa- A possible correla- 
tion between the reddening and L*/Lagn suggested that this dust 
resides on host galaxy scales, beyond the NLR. This suggestion can 
now be tested using Fig. 11, which shows the mean BPT positions 
by Luv/i-bHa- The mean positions of the Luv/^bHa = 40, 10,4, 1 
bins are all similar to each other, in both high luminosity Tls (lower 
row) and low luminosity Tls (upper row). If the dust resides inside 



SF-non robust 69 42.8 0.16 ±10% 37 ±10% 

L matched pure- AGN 207 42.7 0.06 ±10% 29 ±10% 

Table 4. The mean L^ua /LbHa and Luv /LbHa of different BPT-[N II] clas- 
sifications. The SF objects are divided according to whether their nar- 
row line measurements are robust (§2.2.3). In the pure-AGN and compos- 
ite groups, 83% of the objects have robust measurements. The LbHa ( m 
log erg s~'), L,Ha/L>Ha and ^uv/LbHa values are the geometrical means, 
with the uncertainty in the mean noted only if it is > 5%. 



to 0.12 in the matched pure-AGN. Therefore, the L nH a/AjHa ratios 
of Composites are consistent with a roughly equal AGN and host 
contribution to L n H«. In the robust SFs, the host contribution is 
twice the AGN contribution. An intermediate ratio is seen in the 
non-robust SFs. 

A similar effect is expected in Luv/A)Ha> as star formation 
will contribute only to Luv- Indeed, the mean Luv/A>Ha of Com- 
posites and robust SFs is 50% higher than in the respective matched 
group (25% difference in the non robust SFs). Is the observed 
increase in Luv/^bHa consistent with the observed increase in 
L n Ha/LbHa? Star forming galaxies have a mean Luv /i-nHa = 120 
(Kennicutt & Evans 2012). The Composites show an increase of 
0.15 in L n Ha/LbHa, and are thus expected to show an increase 
of 120 x 0.15 = 18 in Luv/^bHa> which is indeed observed (54 
from 36, Table 4). The robust SF group show an increase of 0.3 
in L n Ha/LbHa, and are thus expected to show an increase of 36 in 
Luv/A>Ha> which is 50% larger compared to the observed rise of 
22. However, the difference is probably consistent within the larger 
uncertainties in this group. In the non robust SF group the expected 
rise in Luv/^bHa is 12, versus an observed value of 9, again con- 
sistent with the uncertainties. 

To summarize, the Tl AGN which reside in the Composites 
and SF regions of the BPT diagrams, also show higher Luv /LbHa 
and L n Ha/LbHa ratios, compared to pure-AGN. In addition, the ra- 
tio of the increase in L n H« and in L vv is consistent with Luv /i-nHa 
observed in star forming galaxies. Thus, AGN powered by accre- 
tion onto a massive BH do not produce SF or Composite line ratios, 
and measurements of such line ratios in AGN implies host contam- 
ination. 

Could host contamination also affect line ratios within the 
pure-AGN regime? Could some of the spread in the BPT diagrams, 
also within the pure-AGN regime, be caused by host contamina- 
tion? Figure 11 presents the mean BPT positions of Tls binned 
by Luv/AjHa- We split the Tl sample to LbHa < 10 42 erg s _1 
(upper panels), and LbHa > 10 42 er g s ~' (lower panels). The lu- 
minosity cut is set where the host contribution to Luv starts to be 
significant (Paper I). At LbHa > 10 42 erg s _1 , objects within the 
Luv /LbHa ^ 40 bins have similar mean positions, but the high- 
est bin Luv /LbHa = 80 is shifted towards the Composite region. 
A similar behaviour is observed at LbHa < 10 42 erg s _1 . Objects 
within the Luv /LbHa ^ 30 bins have similar mean positions, but 
the Luv /LbHa = 90 bin is shifted towards the Composite region. 
The highest bin here has Luv/^bHa = 200, and its mean position is 
within the Composite region. Thus, not only that Composite AGN 
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Figure 11. The mean BPT positions by Luv/i-bHa. for low and high luminosity Tls. Mean positions of 0.5 decade Luv/^bHa bins are denoted by triangles, 
with marker size increasing with Luv/^bHa (mean value noted). The error bars indicate the uncertainty in the mean position. Only bins with > 10 objects are 
shown. The positions of individual objects are shown as gray dots. The solid, dashed and dashed-dot classification lines are as in Fig. 2. Top panels Tls with 
log LbHa < 42. The mean position of objects with Luv/ibHa = 90 and 200 are offset towards the composite region, indicating the excess UV originates from 
star formation in the host galaxy. Bottom panels Tls with log L biia > 42. The Luv/AiHa = 40, 10,4, 1 bins have similar mean BPT positions, indicating a 
similar intrinsic ionizing spectrum at these different Luv/^bHa- Therefore, if Z-uv/i-bHa < 40 indicates dust extinction (Paper I), the extincting dust resides on 
scales larger than the NLR. 



the NLR, then the NLR in objects with a low Luv/^bHa is illu- 
minated by a modified ionizing SED, which will shift their mean 
BPT position. The complete lack of a trend in BPT position with 
reddening suggests that the NLR illumination is not modified, and 
therefore the extincting dust resides on scales larger than the NLR. 



9 CONCLUSIONS 

The narrow line ratios of type 2 AGN have been extensively ex- 
plored, in particular based on the SDSS sample. Here we present 
a similar analysis of the Tl sample, the first large (3 175 objects) 
sample of type 1 AGN (Paper I). The Tl sample extends to lumi- 
nosities well below the SDSS quasar sample, and thus in contrast 
with quasars, where the narrow lines are generally difficult to mea- 
sure, here a significant fraction of the objects have strong narrow 
lines (Paper II). This allows reliable analysis of the narrow line ra- 
tios for most objects, as done in type 2 AGN. We find the following: 

(i) The luminosities of all measured narrow lines, Ha, H/3, 
[O III], [N II], [S II], [O I], show a Baldwin relation rela- 
tive to the broad Ha luminosity, L]; ne <>= ^bHcr w i m a ~ 
0.66, 0.67, 0.72, 0.54, 0.53, 0.63, respectively (Paper II, and 
above). 

(ii) About 20% of the Tl AGN have line ratios within the 'Composite' 



and 'SF' regions of the BPT diagrams. These line ratios are not 
powered by accretion onto a massive BH, as these objects also show 
higher L n ua/^bHa an d ^uv/^bHa emission. The excess L n Ha an d 
Lyv is consistent with the ratio expected from SF in the host galaxy, 
and indicates the line emission in these objects is mostly excited by 
SF, rather than by the AGN. 

(iii) The other 80% of the Tl AGN, which reside within the BPT AGN 
region, are offset to lower [N Il]/Ha and [S Il]/Ha luminosity ra- 
tios, compared to type 2 AGN. This offset is a selection effect, as 
Tl AGN selected only from the SDSS galaxy sample, as the type 
2 AGN are, are not offset. The offset is produced by the Tl point 
like objects, selected from the SDSS quasar sample, which extend 
to higher L^gn/^*- The T2 sample is selected against such objects, 
and such objects are offset to lower [N Il]/Ha and [S Il]/Ha. 

(iv) The [N Il]/Ha and [N II]/[S II] ratios increase with host mass, 
which suggest a mass-metallicity relation in AGN hosts, as ob- 
served in quiescent galaxies. In contrast [N Il]/Ha decreases with 
Lagn, but [N II]/[S II] is independent of Lagn, which indicates 
there is no direct LAGN- meta lli c ity relation. 

(v) At a fixed Lb i and L/L E ^, objects with a higher [N Il]/Ha also 
have higher broad Fe II luminosity, suggesting the broad line metal- 
licity is also related to the host mass. This may be an additional 
independent effect to the L/LEdtj-metallicity relation, suggested in 
earlier studies to explain some of the eigenvector 1 relations. 
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Figure 12. The mean BPT-[0 I] positions of the Tl sample, for differ- 
ent a cx . The solid, dashed and dashed-dot classification lines are as in 
Fig. 2. To avoid host contamination of Luv, only Tl objects with L bHa > 
jq42.5 erg s -l use( j (-g ra y dots). 7 ne mean position and scatter of each 
a x bin is marked by a triangle with error bars (mean a ox noted). For com- 
parison, the expected BPT position for ionizing spectra with different slopes 
and different log U are marked (values from Groves et al. 2004). The mean 
BPT positions of the Tl objects do not follow the trend expected if the slope 
of the ionizing spectrum is a ox . Either a ox does not represent the EUV spec- 
tral slope, or other parameters, such as metallicity, ionization parameters 
and density control the scatter. 



(vi) The fraction of AGN which are LINERs increases from ~ at 
V L Edd = 10~ 2 to - 1 at L/L Edd = 1(T 4 . The Tl and T2 samples 
show a similar fraction at a given L/LEdd, indicating the LINER 
phenomena is unrelated to the presence of an observable BLR. 
However, the Tl sample terminates at L/L^dd ~ 10 3 > either due to 
a physical effect or due to selection effects, and thus LINERs con- 
stitute only ~ 4% of the Tl sample, but ~ 50% of the T2 sample. 

(vii) The BPT position is unaffected by the value of Luv /LbHa f° r val- 
ues < 30, which provide a measure of the foreground dust extinc- 
tion (Paper I). This suggests that the ionizing continuum observed 
at the NLR is unaffected by dust extinction, and the dust likely re- 
sides on the host galaxy scale. 

(viii) The BPT position is unaffected by the observed a ox . Models show 
there is a strong dependence of the BPT position on the ionizing 
continuum slope. This suggests that a ox is not a good indicator of 
the EUV ionizing spectrum slope, or that other parameters, such as 
metallicity, density, and ionization parameter, dominate the scatter 
in the BPT plots. 

We thank Dan Maoz and Hagai Netzer for helpful suggestions 
and comments. We thank God for finishing this work. 
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Appendices 

APPENDIX A: MEAN SPECTRA 

In §§3 - 5, we show that the distribution of [N Il]/Ha in the Tl 
sample shifts to lower values with increasing L^ea- To exclude the 
possibility that this trend is an artifact of our deblending algorithm, 
we examine the mean spectra at different L^ea and Av. 

We divide the Tl objects with [O IIl]/Hj3 > 5 to bins of 0.3 
dex in Av and one decade in L^ a . For each bin, we derive the 
host-subtracted mean spectrum, as described in §5.3.2. These mean 
spectra are plotted in Figure Al, with Lbe a increasing from bottom 
to top, and Av increasing from left to right. The thin black line in 
each panel shows the relevant mean spectrum in the vs. velocity 
v plane, centered on Ha. To enhance the contrast between different 
NLR and BLR components, we also plot the same mean spectrum 
centered on H/3 (thick gray line), with the L x of the H/3 region fit 
to the of the Ha region. The fit is performed by a least square 
minimization of two parameters a and b so that 

L\(vna) ~axL A (v H|8 )+Z? (Al) 

for all v in the ranges -3000 < v < -1000 km s _1 , -200 < v < 
450 km s" 1 and 1600 < v < 3000 km s _1 . The v ranges for the 
fit are chosen to avoid the [N II] lines. The best fit coefficient a is 
between 2.7 and 4 in all panels. 

The decrease in [N Il]/Ha with increasing L^ea is clear in the 
Av = 2800, 5200, and 9600 km s columns, confirming the trend 
found on single objects in §§3-5. 



APPENDIX B: BPT BY M BH AND M* 

Figures Bl and B2 show the BPT positions of the Tl sample, di- 
vided by M BH and M«, in the same format as Figures 3-5 above. In 
Fig. B2, only the 91% of Tl objects with a reliable estimate of M t 
(§2.3.1) are shown. 

With decreasing M BH the fraction of Tls classified as Com- 
posites and SFs increases, indicating an increase in the relative 
amount of host contribution to the NLR (§7). Since the SDSS is 
a flux limited sample, Tls with low M BH , and therefore low bulge 
mass, are preferentially selected from disk dominated galaxies (see 
Fig. 16 in Paper I). Disks have a relatively large specific star forma- 
tion rate, which may cause the observed shift in the BPT positions. 

With decreasing M„ an increasing fraction of objects are off- 
set to low [N Il]/Ha values, as found by Groves et al. (2006) on 
a type 2 AGN sample. This trend is consistent with the M t — Z 
relation of quiescent galaxies (Lequeux et al. 1979, and citations 
thereafter). 
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VBalmer (1000 km S ') 



Figure Al. Mean host-subtracted spectra of the [O III]/Hj3 > 5 Tl objects, near Ha and near Hj8, at different LbHa and Av. Each panel shows the mean 
spectrum of Tl objects with the same Av and the same L bHa (number of objects noted). The LbHa increases from bottom to top (mean LbHa in erg s" 1 noted 
in left column), and Av increases from left to right (mean Av in km s _1 noted in top row). The thin black line plots the mean spectrum centered around Ha. 
To enhance the contrast between different NLR and BLR components, we also plot the mean spectrum centered around Hj8, with the L^ adjusted to fit the L x 
of the Ha region (thick gray line). In the Av = 2800, 5200, and 9600 km s~' columns, the [N II]/Ha ratio clearly decreases with increasing LbHa- 



20 Jonathan Stern and Ari Laor 




log [Nil] /Ha log [SII]/Ha log [OI]/Hce 

Figure Bl. As in Fig. 3, for the dependence of BPT position on Mbh- In each row, Tl AGN from a given decade-wide bin in Mbh are plotted (mean Mbh 
noted, in M,..). The frequency of composites and SFs decreases with increasing Mbh, from 32% at log Mbh — 6.3, to 23%, 14% and 6% at log Mbh = 7.1, 7.9 
and 8.8, respectively. 
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Figure B2. As in Fig. 3, for the dependence of BPT position onM t . In each row, Tl AGN from a given half decade wide bin in M, are plotted (meanA/„ noted, 
in Mq). With decreasing M„, an increasing fraction of objects are offset to low [N II]/Ha values, consistent with the M, —Z relation of quiescent galaxies. 



